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Abstract 

Thermoelectric  waste  heat  recovery  is  investigated  for  current  thermoelectric  materials  with  advanced 
heat  exchangers.  Numerical  heat  exchanger  models  integrated  with  models  for  Bi2Te3  thermoelectric 
modules  are  validated  against  experimental  data  from  previous  cross  flow  heat  exchanger  studies  as  well  as 
experiments  using  thermoelectrics  between  counterflow  hot  water  and  cooling  air  flow  channels.  The 
models  are  used  in  optimization  studies  of  thermoelectric  waste  heat  recovery  with  air  cooling  in  a  cross 
flow  heat  exchanger.  Power  losses  from  an  air  fan  and  a  fluid  pump  result  in  an  optimal  configuration  at 
intermediate  cooling  air  and  hot  fluid  flows.  Results  show  that  heat  exchangers  with  Bi2Te3  thermoelectrics 
can  achieve  net  power  densities  over  40  W/l. 
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1.  Introduction 

Historically,  power  generation  from  waste  heat  streams  using  thermoelectric  (TE)  materials  has 
been  overlooked  due  to  low  thermal  efficiencies  (typically  rjth  <  4%).  The  two  limiting  factors  in 
TE  waste  heat  recovery  are  TE  conversion  efficiency  (t/TE  =  electric  work  produced/heat  flux 
through  TE  material)  and,  to  a  lesser  extent,  effective  heat  exchanger  design.  t/th  is  related  to  t/TF 
by  the  following  equation: 

1th  =  »/te  (!) 
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Nomenclature 

Abt  total  external  heat  transfer  area  (m2) 

A0  minimum  free  flow  area  (m2) 

Acore  frontal  area  of  heat  exchanger  (m2) 

Cp  specific  heat  capacity  at  constant  pressure  (J/kg  K) 
da iaj,  (5?min  major  and  minor  dimensions  of  hot  fluid  tubes  (m) 
c/tube  hydraulic  diameter  of  hot  fluid  channel  (m) 

/  Fanning  friction  factor 

h  heat  transfer  coefficient  (W/m2  K) 

/te  current  (A) 

j  Colburn  factor 

Ky E  thermal  conductance  of  thermoelectric  couples  (W/K) 

/fin  height  of  fins  in  air  flow  passage  (m) 

/louv  height  of  louvers  (m) 

/te  height  of  TE  couples  (m) 

m  mass  flowrates  of  air  and  hot  fluid  (kg/s) 

«te  number  of  couples  per  unit  surface  area  of  hot  fluid  tube 
Nu  Nusselt  number 

Pfm  pitch  of  fins  (m) 

/?iouv  pitch  of  louvers  on  fins  (m) 

Aube  spacing  between  tubes  (m) 

P  pressure  in  flows  (Pa) 

Pr  Prandtl  number 

Qte  total  rate  of  heat  dissipated  through  thermoelectrics  (W) 
qsc,  qsh  heat  flux  per  unit  area  of  tube  wall  (W/m2) 

R  electrical  resistance  (Q) 

Re  Reynolds  number 

T  temperatures  of  air  flow,  hot  fluid  flow,  and  thermoelectrics  (K) 
Tsc,  Tsb  cold  side  and  hot  side  temperatures  of  TE  legs  (K) 
tfin  fin  thickness  (cm) 

V  volume  flow  rate  of  air  and  hot  fluid  (m3/s) 

^ore  heat  transfer  volume  of  heat  exchanger  (m3) 

flpump,  flfan  power  demanded  by  fluid  pump  or  cooling  air  fan  (W) 

Ate  total  power  generated  by  TEs  (W) 

wte  power  generated  per  unit  area  of  hot  fluid  tube  (W/m2) 

wTE  width  of  individual  TE  leg  (cm) 

x  axial  distance  along  hot  fluid  tube  flow  axis  (m) 

Z  figure  of  merit  of  thermoelectric  (K  ') 

$he  heat  exchanger  costs 

$te  thermoelectric  costs 
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Greek  symbols 

a  Seebeck  coefficient  (V/K) 

e  heat  transfer  effectiveness 

ss  surface  roughness 

X  thermal  conductivity  of  fluids  or  TE  material  (W/m  K) 

f/TE  conversion  efficiency  of  TE  modules 

rjth  overall  efficiency  of  waste  heat  conversion  to  electricity 

p  density  of  flows  (kg/m3)  or  electrical  resistivity  of  TE  material  (Q  m) 

pc  contact  electrical  resistivity  (Q  m2) 

o  ratio  of  minimum  free  flow  area  to  frontal  area 

9  louver  angle  (°) 

p  dynamic  viscosity  (kg/m  s) 

Subscripts 
a  air  flow 

f  hot  fluid  flow 

TE  TE  materials 


where  e  is  the  fraction  of  waste  heat  passed  through  the  TE  couples.  ;?TE  of  a  thermoelectric  power 
generator  depends  on  the  material’s  non-dimensional  figure  of  merit,  ZT  =  (a^E7"/pTE/.TF)  where 
the  average  temperature  of  the  TE  couples  is  T  =  (Tsh  +  Tsc)/2,  and  the  material  properties  a,  p 
and  X  are  averaged  between  the  specific  properties  of  the  n-  and  p-legs  of  the  TE  as  discussed 
in  previous  references  [1].  Analysis  for  rjTE  in  an  optimized  configuration  [1]  reveals  how  >?TE 
improves  by  increasing  Tsh  and  decreasing  Tsc,  as  indicated  in  Eq.  (2). 


"^+(i  +  -'M(i-t) 


(2) 


where  m'  is  the  ratio  of  the  load  resistance  to  the  total  internal  resistance  of  the  TE  couples.  rjTE 
can  be  readily  shown  to  be  a  monotonically  increasing  function  of  ZTsh.  For  maximum  efficiency, 
m'  =  (1  +  Z(Tsh  +  Tsc)/2)]/2,  and  for  maximum  power  output,  m!  =  1  [1,2].  The  thermal  efficiency 
of  the  waste  heat  recovery  process  can  then  be  evaluated  in  combination  with  heat  transfer 
analysis  to  determine  e. 

Recent  advances  in  materials  and  materials  processing  as  summarized  in  recent  reviews  [3-5] 
have  led  to  higher  ZT  values  and,  thus,  higher  theoretical  f/TE.  For  large  Tsh  —  Tsc  (>500  K),  TE 
modules  combining  various  materials  in  long  segmented  couples  have  exhibited  riTE  near  15%  or 
more  for  7^  ^  850  K  [6,7].  These  multi-component  TE  couples  may  have  applications  in  designs 
where  relatively  high  temperature  heat  sources  are  available  and  where  TE  size  does  not  greatly 
impact  heat  exchanger  effectiveness.  For  applications  where  Tsc  is  well  below  room  temperature, 
ZT  up  to  2.4  [8]  have  been  discussed  for  Tsh  near  room  temperature.  Very  high  ZTs  have  also  been 
discussed  for  proposed  quantum  well  TE  devices  [9]  such  that  j?te  may  approach  20%  for  some 
applications. 
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For  higher  temperature  applications,  new  TE  materials,  such  as  Zn4Sb3  based  semiconductors 
[5]  and  skutterudites  [3,5,6,10],  offer  improved  f/TE,  much  higher  than  the  5%  of  Bi2Te3,  but  to 
achieve  such  performance  with  these  new  materials  TsU  must  be  >600  K.  For  waste  heat  recovery 
applications,  such  high  temperatures  require  direct  heat  transfer  from  gas  phase  exhaust  streams 
or  the  use  of  expensive  high  temperature  liquids,  such  as  liquid  metals.  In  general,  direct  heat 
transfer  from  gas  phase  exhaust  flows  will  result  in  either  a  low  e  or  a  large  power  penalty  in  terms 
of  pressure  drop  across  the  heat  exchanger.  To  operate  effectively  for  waste  heat  recovery 
applications  where  inexpensive  fluids,  such  as  water  or  water/glycol  mixtures,  are  used  to  extract 
the  heat  from  exhaust  flows  (as  in  automotive  radiators  or  most  cooling  tower  applications),  TE 
materials  must  operate  effectively  at  350  K  <  Tsh  <  400  K.  With  presently  available  materials, 
Bi2Te3  based  semiconductors  offer  the  best  performance  for  this  temperature  range  with  ZT 
approaching  1.1  [5].  For  Tsh  ~  120  °C  and  Tsc  «  20  °C,  this  ZT  gives  >?TF  up  to  5%  from  Eq.  (2). 
The  relatively  affordable  Bi2Te3  TEs  are  still  one  of  the  more  attractive  options  for  waste  heat 
recovery  and  are  the  focus  of  the  current  study. 

Previous  experiments  have  investigated  waste  heat  recovery  using  Bi2Te3  TEs  and  have  shown 
overall  waste  heat  recovery  efficiencies  ( rith )  <  4%  [11,12],  Examination  of  the  results  reveals  that 
the  heat  exchanger  design  on  both  the  hot  and  cold  sides  of  the  TEs  significantly  impacts  rjth.  The 
current  study  addresses  issues  related  to  optimizing  TE  waste  heat  recovery  by  integrating  efficient 
cross  flow  heat  exchangers  with  thermoelectric  modules  for  conversion  of  waste  heat  to  electricity. 
Numerical  heat  exchanger  models  are  integrated  with  models  for  Bi2Te3  TE  modules,  which  are 
positioned  between  the  hot  and  cold  flow  streams  as  shown  in  the  schematic  of  a  heat  exchanger 
sub-section  in  Fig.  1.  Fig.  1  indicates  fins  on  the  air  side  passages,  which  are  often  louvered  (not 
shown  in  the  illustration)  to  enhance  heat  transfer,  since  air  side  cooling  typically  presents  the 
highest  resistance  to  heat  transfer. 


cooling  air  flow 


Fig.  1.  Schematic  showing  layout  of  sub-section  of  thermoelectric  heat  exchanger. 
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The  numerical  model  incorporates  the  systems  level  power  requirements  associated  with  both 
pumping  the  heat  source  (liquid)  and  blowing  the  cooling  air  through  the  exchanger.  The 
numerical  model  is  validated  in  two  parts.  The  heat  exchanger  model  is  validated  against  mea¬ 
sured  performance  of  advanced  cross  flow  heat  exchangers  without  thermoelectrics,  and  the  part 
of  the  model  with  integrated  thermoelectric  modules  is  validated  against  simple  experiments  using 
thermoelectrics  between  counterflow  hot  water  and  cooling  air  flow  channels.  The  full  model 
is  used  in  a  systems  level  optimization  to  determine  the  optimal  performance/design  of  thermo¬ 
electric  power  generators  for  waste  heat  recovery  from  hot  water/glycol  mixture  streams.  Opti¬ 
mization  of  a  thermoelectric  power  generation  module  has  been  previously  presented  [13], 
but  these  optimization  studies  did  not  address  system  accessories  and  constraints  associated 
with  the  external  heat  transfer.  The  previous  studies  indicated,  for  a  given  design,  longer  TE 
couples  tend  to  produce  higher  rjTE  but  at  the  expense  of  lower  e  and,  in  general,  lower  tjih.  Other 
studies  have  incorporated  power  penalties  associated  with  the  hot  fluid  flow  pressure  drops 
[14,15],  and  these  studies  identified  some  of  the  trade  offs  between  TE  couple  geometry  and 
efficient  heat  exchange.  This  new  study  focuses  on  simultaneously  optimizing  both  the  heat  ex¬ 
changer  geometry  as  well  as  the  TE  geometry  for  optimal  net  waste  heat  recovery  performance  by 
including  the  penalties  associated  with  both  the  hot  and  cooling  flow  pressure  drops.  This  study 
also  provides  a  basis  for  further  analysis  as  new  TE  materials  and  segmented  modules  become 
available. 


2.  Model  description 

The  basic  configuration  for  the  optimization  study  follows  the  general  design  suggested  by  the 
heat  exchanger  subsection  illustrated  in  Fig.  1 .  The  TE  modules  are  located  along  the  perimeter 
on  the  major  axis  of  the  hot  flow  passages  as  in  Fig.  1.  The  TE  modules  are  electrically  isolated  by 
thin  passivating  oxide  layers  on  both  the  outer  surface  of  the  metal  tubes,  which  hold  the  hot  fluid, 
and  the  inner  surface  of  the  outer  metal  wall,  which  supports  the  fins  for  air  side  cooling.  Inte¬ 
gration  of  TEs  in  the  heat  exchanger  design  must  not  cause  a  significant  drop  in  e,  and  this  re¬ 
quires  that  the  length  of  the  TE  legs  remain  relatively  small.  By  fixing  this  basic  configuration  and 
varying  the  heat  exchanger  and  TE  geometrical  parameters  as  well  as  the  heating  and  cooling  flow 
rates,  significant  insight  is  gained  into  the  trade  offs  and  principles  for  optimizing  TE  waste  heat 
recovery  on  a  systems  level. 

A  numerical  model  has  been  developed  for  a  single  tube  with  hot  fluid  flow,  surrounding  TE 
modules  and  external  cooling  air  flow  over  extended  heat  transfer  surfaces  (i.e.  this  case,  louvered 
fins).  The  model  domain  includes  half  the  width  of  the  external  cross  flow  passages  on  both  sides 
of  the  tube.  Thus,  under  the  assumption  that  all  hot  fluid  flow  tubes  are  performing  equally,  the 
mid-plane  boundary  in  the  cross  flow  passage  between  adjacent  tubes  acts  as  an  adiabatic  surface, 
and  a  single  tube  model  can  capture  the  performance  of  the  entire  thermoelectric  heat  exchanger. 
For  optimization  studies,  the  dimensions  of  the  tubes  and  the  extended  surface  features  in  the 
cross  flow  cooling  passages  are  free  to  vary  with  the  number  of  tubes  in  a  fixed  volume  to 
determine  optimal  net  power  densities.  The  model  is  set  up  with  the  option  to  exclude  the  ther¬ 
moelectric  modules  surrounding  the  heat  exchanger  tube  in  order  to  validate  the  cross  flow  heat 
transfer  correlations  against  experiments  of  similarly  designed  heat  exchangers.  For  this  study,  the 
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heat  exchanger  design  is  constrained  to  a  single  row  of  tubes,  rather  than  multiple  rows  of  inline 
or  staggered  tubes. 


2.1.  Heat  exchanger  modeling 


In  order  to  facilitate  parametric  design  optimization,  the  heat  exchanger  model  implements 
empirical  correlations  for  the  hot  side  and  cold  side  convective  heat  transfers.  For  the  hot  fluid 
flow  through  the  inner  tube,  the  Nusselt  number,  Nut  =  h,‘{"‘hc,  uses  correlations  for  fully  developed 
pipe  flow,  since  the  ratio  of  /tube  to  hydraulic  diameter  dtube  typically  ranges  from  40  to  100.  For 
the  laminar  flow  regime,  Nu{,  assuming  constant  wall  temperature,  varies  from  2.9  to  7.5  with 
variations  in  the  aspect  ratio  of  the  tubes  (dmaj/dmin)  [16].  For  transition  and  turbulent  flow 
conditions,  Gnielinski’s  correlation  for  Nusselt  number  is  used,  where  the  property  values  are 
evaluated  at  the  local  fluid  temperature  T{  [16]. 


(ft/2)(Rei  -  1000)iVf 
1  +  2.7(/f/2)1/2(FVf/3  —  1) 


(3) 


The  Fanning  friction  factor  fi  for  the  tube  flow  is  calculated  for  all  flow  regimes  using  the  fol¬ 
lowing  formula  developed  by  Bhatti  and  Shah  [17] 


ft 


=  A  + 


B 

Re^ 


(4) 


where  for  laminar  flow  (Ref  <  2100),  A  =  0,  B  =  16.0  and  m  =  1.0,  for  transition  flow  (2100  < 
Ref  ^4000),  A  =  0.0054,  B  =  2.3  x  10  8  and  m  =  —0.6667,  and  for  turbulent  flow  (Ref  >  4000), 
A  =  0.00128,  B  =  0.1143  and  m  =  3.215.  Using  this  Fanning  friction  factor,  the  fluid  side  pressure 
are  drop  and  associated  pumping  work  are  calculated. 

For  the  external  cooling  air  flow,  the  heat  transfer  coefficient  over  the  extended  surface  uses  a 
previously  developed  correlation  [18] 


K 


jmaCp,a 

Ao#? 


(5) 


where  the  Colburn  factor,  j,  is  taken  from  a  generalized  correlation  for  a  corrugated  louver  fin 
geometry  for  cross  flow  heat  exchangers  [19]. 


where  Rea  is  calculated  with  the  velocity  at  A0  and  a  length  scale  of  p\ ouv.  This  correlation  was 
validated  from  a  compilation  of  many  studies  for  a  range  of  louver  geometries  and  shows  good 
agreement  with  experiments  (±25%  for  98%  of  the  samples  compared)  over  a  broad  range  of  flow 
conditions  with  Rea  ranging  from  100  to  3000.  Although  other  correlations  exist  for  flow  over 
louvered  fins  [20-23],  the  above  correlation  was  chosen  for  the  extent  of  its  validation  and  its  more 
general  nature. 

The  cooling  air  flow  pressure  drop,  which  determines  the  required  cooling  fan  power,  is  cal¬ 
culated  by  using  [24] 
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APa  =  - 


\  Pa, out  /  ^0  Pa, mean 


"  2^0Pa,i 
where  /a  is  based  on  Red  [25] 

/a  =  0.805i?e~°514  (  (  — 

VPlouv  /  VPlouv  , 


\Plouv  / 


(7) 

(8) 


The  friction  factor  correlation  is  valid  for  Rea  =  100  —  700  and  ratios  of  louver  fin  surface  area  to 
outer  tube  surface  area  ranging  from  7  to  12.  Other  friction  factor  correlations  for  different 
louvered  fin  geometries  exist  [23],  but  the  selected  correlation  provides  consistency  with  the  heat 
transfer  correlation.  A  comprehensive  analysis  comparing  different  louver  fin  designs  has  not  been 
completed  for  this  study. 


2.2.  Thermoelectric  modeling 

As  stated  in  the  introduction,  thermoelectric  modules  of  Bi2Te3  are  well  suited  for  applications 
for  low  temperature  waste  heat  recovery.  The  TE  properties  for  both  p-type  and  n-type  Bi2Te3 
semiconductors  are  calculated  by  using  information  provided  by  Melcor,  Inc.  [26],  A  simple 
energy  balance  around  the  TE  modules  yields  the  following  equations  for  the  hot  side  heat  flux 
into  the  TEs,  the  cold  side  heat  flux  out  of  the  TEs  and  the  power  generated  per  unit  surface  area 
of  the  hot  fluid  tube: 

</sh  =  (KTE(Tsh  -  Tsc)  +  <xTshI  -  1 72i?TE  -  I2Rc^j  *  nTE  (9) 

<7sc  =  ^te(7*  -  Tsc)  +  rj.TJ  +  I/2i?TE  +  l2RL)j  *  „TE  (10) 

h’TE  =  («(Tsh  —  Tjc)/  —  /2(7?te  +  2 Rc))  *  «te  (11) 

The  material  properties  for  the  thermoelectrics  in  Eqs.  (9)— ( 11)  vary  along  the  length  of  the  hot 
fluid  tube  with  changes  in  temperature.  The  electrical  contact  resistance  Rc  between  the  p  and  n 
couples  was  taken  into  account  as  in  previous  references  [2,27], 

Thermal  contact  resistance  arises  between  the  TE  couple  interconnects  and  the  oxide  layers  on 

the  heat  exchanger  material.  It  is  taken  into  account  by  simulating  a  partial  air  gap,  which  in¬ 

creases  the  resistance  to  heat  transfer  and  causes  a  small  drop  in  Tsh  between  the  two  adjoining 
surfaces.  These  air  gaps  represent  the  surface  irregularities  in  the  layers,  which  increase  the 
resistance  to  heat  transfer  and  thereby  reduce  the  temperature  difference  across  the  TE  couples. 

The  current  flow  through  the  thermoelectrics  for  a  given  side  wall  of  a  heat  exchanger  tube 
is  determined  by  the  number  of  couples  aligned  in  series. 

J  E”=l  («(rsh  -  TSC)»TE  *  ^maj  Ax)  . 

TE  ^ad  +  EE  ((/?n:  +  2Rc)nn:  *  d^Ax),  1  J 

where  n  =  /tube/At.  The  numerator  in  Eq.  (12)  is  the  sum  of  the  voltages  generated  in  each  z'th 
segment  of  the  axial  discretization,  and  the  summation  in  the  denominator  represents  the  internal 
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resistance  due  to  the  couples  and  contact  resistances  in  each  segment.  Eq.  (12)  applies  when  the 
thermoelectric  modules  are  wired  in  series  along  each  tube  but  in  parallel  for  separate  tubes.  The 
power  per  tube  is  simply  /2i?load. 

2.3.  Numerical  approach 

The  numerical  model  is  implemented  in  MATLAB.  The  primary  discretization  is  along  the 
axial  direction  of  the  hot  fluid  flow,  which  is  captured  with  a  first  order  upwind  differencing 
scheme.  The  transverse  heat  transfer  out  of  the  hot  flow  is  modeled  using  the  above  heat  transfer 
correlations  and  TE  analysis.  Thermal  conduction,  both  across  and  along  the  walls  of  the  tubes,  is 
incorporated  using  central  differencing  discretization  for  the  gradients.  For  the  air  side  cross  flow, 
the  log  mean  temperature  is  solved  for  at  each  axial  location  and  is  used  to  determine  the  heat 
transfer  to  the  cooling  air  flow. 

The  cells  are  set  at  a  constant  length  Ax,  specified  as  a  multiple  of  the  distance  between  adjacent 
rows  of  thermoelectric  couples.  For  the  scenarios  in  this  paper,  Ax  ranges  from  1  to  2  cm.  In  light 
of  the  mild  temperature  gradients  across  the  heat  exchanger,  this  level  of  discretization  is  deter¬ 
mined  to  be  more  than  adequate  for  accurately  modeling  the  variation  in  heat  fluxes  and  ther¬ 
moelectric  power  across  the  tube  length.  The  discretized  energy  conservation  equations,  as 
discussed  in  detail  in  previous  references  [28,29],  lead  to  equations  at  all  axial  distances  along  the 
tube  flow  for  T{,  Tsh,  Tsc,  T  in  the  inner  and  outer  tube  walls  and  the  log  mean  of  T.d.  The  model 
solves  these  non-linear  governing  equations  simultaneously  for  the  steady  state  temperatures  via 
an  iterative  Newton-Raphson  method  using  an  analytical  Jacobian. 

Along  with  the  temperatures,  the  current  flow  through  the  thermoelectric  couples,  7,  must  be 
determined.  Since  the  couples  for  each  tube  are  wired  in  series,  as  indicated  by  Eq.  (13),  I  is 
constant  for  all  axial  positions.  Numerically,  the  global  nature  of  I  causes  the  Jacobian  of  the 
governing  equations  to  lose  its  block  tri-diagonal  structure,  which  slows  the  Newton-Raphson’s 
iterations,  although  not  prohibitively. 


3.  Experimental  validation 

3.1.  Thermoelectric  heat  transfer  validation  experiment 

An  experiment  was  set  up  to  validate  the  heat  transfer  model  through  the  thermoelectrics  as 
detailed  above.  A  block  diagram  of  the  experiment  is  shown  in  Fig.  2.  A  single  aluminum  tube  (40 
mm  x  8  mm  x  0.5  mm),  which  flowed  hot  water,  was  lined  on  each  side  with  three  TE  modules  (0.4 
cm  thick)  that  included  thin  A1203  plates  around  the  TE  couples.  On  the  cold  side  of  the  TE 
modules,  aluminum  tubes  for  the  cooling  air  flows  were  pressed  against  the  TE  ceramic  insulating 
plate  (as  shown  in  Fig.  2).  Silver  based  thermal  grease  (manufactured  by  Arctic  Silver  Inc.)  was 
placed  between  all  of  the  ceramic  plate/flow  tube  interfaces  in  order  to  minimize  the  thermal 
contact  resistance. 

Hot  water  through  the  center  tube  was  pumped  around  a  short  closed  loop  system  that  included 
a  reservoir  with  a  temperature  controlled  immersion  heater  and  a  flow  meter  as  indicated  in  Fig.  2. 
To  ensure  even  distribution  of  the  water  flow  through  the  tube,  flow  straighteners  were  used  at  the 
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Fig.  2.  Schematic  of  counterflow  TE  heat  exchanger  experiments  used  for  TE  model  validation. 


exit  to  backpressure  the  water  flow  and  force  it  to  wet  the  entire  tube  surface.  Electronic  mass  flow 
controllers  controlled  the  high  air  flows  through  the  cooling  tubes.  The  air  flows  in  the  coun¬ 
terflow  heat  exchange  configuration  provided  high  cooling  and  permitted  the  use  of  a  well 
established  heat  transfer  correlation  in  modeling  the  experiments.  This  allowed  the  modeling  of 
this  validation  experiment  to  focus  on  the  models  for  the  TE  module  effective  conductivity  and 
contact  resistances.  Temperature  measurements  were  made  using  K-type  thermocouples  at  the 
inlet  and  outlet  of  each  of  the  tubes.  In  addition,  thermocouples  were  placed  at  the  center  of  each 
of  the  ceramic  surfaces  of  the  TE  modules. 

Heat  leakage  from  the  outer  tube  walls  forced  the  experimental  model  to  incorporate  heat 
losses  from  the  minor  dimension  on  the  hot  fluid  flow  tube  and  from  the  outer  walls  on  the  cooling 
air  flow  tubes.  The  heat  losses  were  determined  through  fits  used  to  match  the  changes  in  tem¬ 
perature  and  power  output  as  a  function  of  load  resistance.  Although  the  heat  losses  were  kept 
small  with  alumina-silicate  fiber  insulation  around  the  rig,  modeling  these  losses  was  necessary  to 
capture  the  observed  (Tsh  -  Tsc).  Computer  simulations  from  the  revised  model  were  compared 
with  experimental  results  for  the  configuration  shown  in  Fig.  2.  A  least  squares  optimization  was 
used  to  determine  the  heat  transfer  coefficients  for  the  external  heat  losses  and  to  fit  the  thermal 
contact  resistances  between  the  ceramic  layers  around  the  TE  couples  and  the  metallic  tubes.  In 
addition,  the  electrical  contact  resistivity  pc  was  determined  by  a  least  squares  fit  of  the  experi¬ 
mental  data  to  the  model  predictions.  A  value  of  5  x  10  12  Q  m2  was  found  to  give  good  agreement 
with  the  contact  resistances  that  have  been  measured  and  reported  in  the  literature  [30,31].  Over 
35  experimental  cases  were  performed  for  a  range  of  hot  fluid  and  cooling  air  flows  and  for  a  wide 
range  of  TE  power  conditions,  such  that  the  results  could  be  used  to  evaluate  the  TE  heat  transfer 
model. 

The  thermal  contact  resistances  were  evaluated  in  terms  of  an  equivalent  air  gap  thickness 
between  the  ceramic  plates  of  the  TE  module  and  the  flow  tubes  on  the  hot  and  cold  sides  of  the 
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TE.  Thicknesses  for  the  equivalent  air  gaps  were  calculated  to  be  only  a  few  hundredths  of  a 
millimeter,  and  the  effective  thermal  contact  resistances  for  the  experimental  conditions  were 
found  to  be  >■*  1 400  W/m2K,  which  also  corresponds  well  with  values  reported  in  the  literature 
[32]. 

Results  from  the  validation  experiments  are  summarized  in  Figs.  3  and  4.  Fig.  3  plots  current 
versus  power  where  the  curves  were  generated  by  varying  the  external  load  resistance,  f?ioad,  at 


0  0.05  0.1  0.15  0.2 

Current  (A) 

Fig.  3.  Experimental  data  and  model  results  for  TE  power  versus  current  for  three  different  sets  of  conditions.  Curve 
was  created  by  taking  measurements  at  various  load  resistances. 


Fig.  4.  Temperature  along  tube  for  the  hot  fluid,  cold  air  and  on  the  hot  and  cold  ceramic  surfaces  of  the  thermo¬ 
electrics  (Vf  =  1.52  1pm,  K,  =  120  slpm). 
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different  fluid  and  air  flows.  Previous  studies  have  reported  that  the  maximum  power  of  a  TE 
module  is  achieved  when  i?ioad  equals  the  effective  internal  resistance  of  the  module  [1,2,27],  The 
calculated  total  internal  resistance  for  the  three  curves  in  Fig.  3  are  33.6  ±  0.4  £1  and  the  peak 
of  the  simulated  power  curves  occurs  where  7?ioad  equals  the  internal  resistance. 

Fig.  4  compares  the  measured  temperatures  with  those  predicted  by  computer  simulation. 
Agreement  is  within  ±2  K  except  for  the  cold  side  temperature  of  the  TE  ceramic.  The  experi¬ 
mental  cold  side  temperature  is  lower  than  the  simulated  temperatures  most  likely  because  the 
experimental  measurement  was  taken  at  the  interface  of  the  TE  ceramic  with  the  surrounding 
metal  tube  and,  thus,  took  an  average  of  the  temperatures  on  the  inside  of  the  cooling  air  tube  and 
the  outside  of  the  TE  module  ceramic  plate. 

A  sensitivity  study  on  the  computational  model  results  was  implemented  to  explore  sources  of 
variation  between  the  experimental  and  model  results.  The  dimensionless  sensitivity  coefficients 
are  highest  for  the  initial  cooling  air  and  hot  fluid  temperatures.  A  1%  change  in  initial  fluid 
temperature  can  cause  as  much  as  a  15%  change  in  power  output  for  some  conditions.  This  large 
sensitivity  of  power  to  fluid  temperature  made  it  critical  to  use  actual  experimentally  measured 
conditions  when  validating  the  model  and  further  suggests  why  the  simulation  curves  for  power  in 
Fig.  3  deviate  from  the  experimental  data  for  some  cases. 

The  relative  magnitudes  of  the  resistances  to  heat  transfer  between  the  hot  fluid  and  the  cooling 
air  in  the  experimental  configuration  can  be  inferred  from  the  temperature  plots  in  Fig.  4.  Even 
with  the  high  air  velocities  in  the  experiments,  the  cold  air  side  heat  transfer  is  not  optimal,  and 
40%  of  the  overall  temperature  drop  occurs  between  the  thermoelectric  module  and  the  cooling 
fluid,  as  suggested  in  Fig.  4.  Nonetheless,  for  the  experiments,  approximately  40%  of  the  total 
temperature  drop  between  7}  and  T.d  does  occur  across  the  TE  module,  which  indicates  that  these 
experiments  provided  a  useful  configuration  for  validating  the  heat  transfer  model  for  the  TEs.  An 
optimal  heat  exchanger  for  waste  heat  recovery  would  maximize  this  temperature  difference  across 
the  TE  couples  by  having  minimal  heat  transfer  resistance  and,  thus,  temperature  drop  at  both  the 
hot  and  cold  side  fluid-wall  interfaces. 

3.2.  Cross  flow  heat  transfer  validation 

The  cross  flow  correlations  for  the  air  side  heat  transfer  were  implemented  and  validated  for 
heat  exchangers  without  thermoelectrics  in  a  similar  manner  as  previously  reported  [19].  Typical 
cross  flow  heat  exchanger  specifications  were  used  as  a  baseline,  both  to  validate  the  model  and  to 
ascertain  what  were  reasonable  limitations  on  dimensions  for  a  design  optimization.  For  this 
study,  the  heat  exchanger  dimensions  were  similar  to  those  of  a  stock  cross  flow  heat  radiator 
(manufactured  by  Harrison  Thermal  Systems  for  a  GM  Chevrolet  Suburban). 

The  cross  flow  heat  exchanger  model  was  validated  against  actual  radiator  performance  data 
without  thermoelectric  modules,  and  some  of  the  results  are  shown  in  Fig.  5a  and  b.  The  model 
predicts  the  axial  drop  in  coolant  temperature,  A Tf,  to  within  4%  for  the  entire  range  of  Va  and  V{ 
as  indicated  in  Fig.  5a.  The  model  also  captures  the  experimental  variation  of  the  pressure  drops 
with  flow  rates,  both  A Pa  with  Va  and  A Pf  with  F{.  The  variation  of  APa  with  Va  is  plotted  in  Fig. 
5b,  and  the  results  indicate  that  the  model  under  predicts  A Pa  by  no  more  than  25%  at  the  highest 
air  flow  rates.  APf  as  a  function  of  V{  is  captured  to  within  5%  over  the  entire  range  of  flow  rates 
and  is,  therefore,  not  shown  in  Fig.  5.  The  model  validation  with  the  radiator  data  provides 
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(a)  Va  (irv'/s) 


(b)  Va  (m3/s) 

Fig.  5.  Results  from  radiator  validation  study  based  on  data  from  Harrison  thermal  systems.  Model  results  indicated 
by  lines  and  experimental  data  indicated  by  points,  (a)  Coolant  temperature  drop  across  radiator  as  a  function  of  air 
flow  over  radiator  tubes  for  a  range  of  coolant  fluid  flows,  (b)  Air  side  pressure  drop  in  cross  flow  as  a  function  of  air 
flows  for  a  range  of  coolant  fluid  flows. 

confidence  that  the  hot  and  cold  flow  heat  transfers  and  pressure  drops  are  captured  well  enough 
that  the  power,  fVTE,  produced  by  the  TE  couples  placed  between  the  hot  and  cold  side  walls  can 
be  calculated  to  a  desirable  degree  of  accuracy.  With  the  validated  heat  transfer  correlations  and 
the  TE  model  parameters,  the  model  was  established  for  optimizing  a  cross  flow  heat  exchanger 
for  TE  waste  heat  recovery. 


4.  Optimization  study 

With  the  cross  flow  heat  exchanger  model  integrated  with  the  models  for  the  TE  couples,  an 
algorithm  was  implemented  to  optimize  the  cross  flow  heat  exchanger  configurations  for  waste 
heat  recovery.  The  objective  function  for  the  optimization  is  the  net  power  (Wnet)  divided  by  the 
heat  exchanger  cost.  Wnet  is  calculated  as  the  total  TE  power  output  from  the  heat  exchanger 
minus  the  expected  power  to  operate  an  electric  fan  for  the  cooling  air  flow  and  an  electric  pump 
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for  the  hot  fluid  (a  water/glycol  mixture).  The  cost  includes  the  thermoelectric  and  aluminum  cost 
of  the  heat  exchanger.  The  aluminum  cost  is  derived  from  cost  rates  for  commercially  available 
radiators,  and  the  cost  of  the  thermoelectrics  is  derived  from  the  bulk  cost  rates  for  commercially 
available  TE  materials.  The  TE  cost  is  divided  into  manufacturing  cost  per  couple,  material  cost 
per  volume  of  TE  material  and  assembly  cost  per  TE  module.  These  costs  are  best  estimates,  and 
the  optimal  solution  will  be  influenced  to  some  degree  by  changes  in  the  cost  functions.  Previous 
studies  looked  at  optimizing  costs  for  integrated  thermoelectric  heat  exchangers  based  on  varying 
fuel  costs  [13,33],  but  unlike  the  current  study,  the  earlier  optimizations  focused  on  the  effects 
of  TE  area  and  leg  length  without  looking  at  the  total  system  cost  and  net  energy  output. 

The  objective  function  for  the  current  optimization,  F(X),  where  X  is  the  vector  of  optimization 
variables,  is  shown  here: 


_  ffnet  _  WTE  -  p  -  Iffan 

($te  +  She)  ($te  +  She) 

(13) 

X  =  [r/maj,  dmm .  /TE,  M-’TE,  Am,  Aube:  K,  V{]T 

(14) 

_  In  order  to  ensure  that  reasonable  solutions  are  found,  bounds  are  placed  upon  the  variables  in 
X.  The  upper  and  lower  bounds  for  all  variables  in  X  are  set  as  discussed  in  detail  elsewhere  [29] 
and,  in  general,  follow  standard  cross  flow  heat  exchanger  geometries  and  expected  manufac¬ 
turing  limitations.  An  additional  constraint  is  set  to  ensure  dmai/dmm  >  5.0  in  order  to  avoid 
having  substantial  amounts  of  heat  loss  through  the  minor  dimension  of  the  tube  where  no  TE 
couples  are  located.  A  final  constraint  implemented  to  ensure  that  the  optimization  scheme  does 
not  converge  to  a  low  cost  heat  exchanger  with  minimal  power  output  by  forcing  the  heat  ex¬ 
changer  to  have  a  Wnet  ^  1.0  kW. 

An  optimized  solution  for  the  cross  flow  heat  exchanger  design  with  a  frontal  area  of  1  mx  1  m 
is  listed  in  Table  1.  The  net  power  output  of  the  optimal  solution  is  the  minimum  of  1  kW,  with  a 
current  of  0.71  A  per  tube  and  a  voltage  of  70.8  V  for  each  tube.  Wiring  each  tube  in  parallel  keeps 
the  voltage  <100  V.  Cost  increases  outweighed  increases  in  Wne t  above  1  kW,  thus  providing  lower 
power/cost  ratios. 

To  assess  the  importance  of  different  parameters  for  optimizing  the  heat  exchanger,  a  para¬ 
metric  study  was  performed  on  the  design  using  the  optimal  solution  as  the  baseline  case.  The 


Table  1 

Optimal  design  for  remote  heat  exchanger  based  on  power/cost  ratio  (kW/$  10,000)  objective  function  and  constant  load 
impedance  of  50  ft 


Design  variables 

Performance  outputs 

dmaj 

2.23  cm 

%e,/($  HE  +  $TE) 

1.09  kW/$  10,000 

drain 

0.45  cm 

Wnet 

1.00  kW 

Ptm 

0.13  cm 

Wte 

1.34  kW 

Aube 

3.74  cm 

Qte 

63.08  kW 

% 

0.17  cm 

Ote 

2.12% 

Wte 

0.20  cm 

Ite 

0.71  A 

K 

2.12  m3/s 

^Internal 

27.4  n 

V{ 

294.1  1pm 
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effects  of  the  variation  of  TE  leg  length  on  net  power  output  and  the  power/cost  ratio  function  F 
are  shown  in  Fig.  6a.  The  peak  at  the  optimal  design  point  occurs  because  of  a  trade  off  between 
riTE  and  e.  Although  higher  t]TE  is  achieved  with  longer  leg  lengths  [13],  the  longer  leg  lengths 
create  more  restriction  to  the  air  flow  between  the  heat  exchanger  tubes  and,  thus,  greatly  reduce 
e.  Previous  studies  of  single  tubes  have  shown  that  when  the  electrical  contact  resistances  are 
negligibly  small,  the  power  output  increases  with  TE  leg  lengths  being  as  short  as  manufacturing 
will  allow.  However,  these  power  gains  come  at  the  expense  of  reduced  rjTE.  On  the  other  hand,  as 
shorter  leg  lengths  reduce  the  electrical  and  thermal  resistances,  the  contact  resistances  begin  to 
become  significant  and  negatively  impact  both  rjTE  and  e  [10,13,32],  The  results  in  Fig.  6a  suggest 
that  for  lTE  <0.15  cm,  the  contact  resistances  begin  to  decrease  the  net  power. 

With  Af  being  nearly  30  times  higher  than  Aa,  ha  is  ~3%  of  ht  .  Thus,  increases  in  heat  transfer 
surface  area  on  the  cold  side  are  needed  to  avoid  the  cold  side  surface  to  air  heat  transfer 
becoming  the  dominant  resistance  to  heat  flow.  The  need  for  increased  air  side  surface  area  is 
reflected  in  the  values  for  optimal  dmai  and  «tube  shown  in  Fig.  6b  and  c.  Increases  in  both  of  these 
parameters  greatly  enhance  the  cold  side  heat  transfer,  but  maxima  arise  in  Wnet  because  A Pa  and 
the  associated  Wian  increase  with  c/ma|  and  «tUbe  due  to  increased  wall  friction  and  decreased  air  flow 
area,  respectively.  A  similar  maximum  in  fVnet  is  found  for  pEm  as  it  decreases  because  the  increased 
number  of  fins  improves  the  air  side  heat  transfer  but  at  the  expense  of  reduced  flow  area  and 
increased  friction  for  the  air  flow.  However,  for  the  optimal  solution  in  Table  1,  pEm  approaches 
the  manufacturing  constraints  (~0. 1  cm),  which  indicates  the  significance  of  increasing  the  heat 
transfer  area  on  the  air  side. 

Fig.  6b  shows  how  the  power/cost  ratio  and  Wne t  vary  differently  with  c/maj.  The  highest  power/ 
cost  ratios  occur  only  when  Wnet  <  1  kW,  and  thus,  the  optimal  point  shown  in  Table  1  is  not  at 
the  highest  possible  value  of  F  due  to  the  power  constraint  of  the  objective  function.  As  dmai 
increases,  the  cost  of  the  exchanger  increases  faster  than  Wnet  because  the  addition  of  more  TE 
couples  (which  is  the  major  cost  factor)  does  not  produce  proportionally  as  much  increase  in 
power. 

The  tube  spacing,  pt ube,  also  has  a  major  impact  on  the  TE  cost.  Therefore,  the  curves  for 
power/cost  ratio  and  Wnet  versus  ptube,  shown  in  Fig.  6c,  do  not  follow  the  same  trends.  Again, 
higher  than  optimal  power/cost  ratios  can  be  achieved,  but  these  values  are  only  found  when  the 
net  power  output  is  below  the  minimum  Wnet  constraint.  The  trade  offs  between  increased  APa  with 
the  associated  increase  in  W{aD  and  increased  fVTE  by  decreasing  pluhe  cause  the  observed  peak  in 
Wne t  with  ptube  in  Fig.  6c.  If  ptube  is  too  high,  the  reduced  number  of  thermoelectrics  results  in 
inadequate  electric  power  to  outweigh  the  losses  for  driving  the  flows. 

For  the  optimization  studies  above,  i?ioad  was  kept  constant  at  50  Q  per  tube.  As  can  be  seen 
from  Table  1,  the  TE  internal  resistance  per  bank  of  modules  on  each  tube  side  is  ~27  Q,  which  is 
approximately  equal  to  Rloai/2.  As  discussed  earlier,  the  theoretical  maximum  power  output  of  a 
thermoelectric  generator  system  occurs  when  the  TE  internal  resistance  equals  R\oa<i-  The  difference 
in  internal  resistance  and  R]0.d(i  at  the  optimal  design  point  arises  because  of  the  penalties  included 
for  costs  and  accessory  power  demand. 

Fig.7  a  shows  how  WTE,  Wpmnp,  and  Wnet  vary  with  V{  for  the  optimal  design  geometry  pre¬ 
sented  in  Table  1.  As  if  increases,  WTE  increases,  but  above  300  1pm,  the  increases  in  WTE  are 
outweighed  by  the  increases  in  Wpump,  and  thus,  there  is  a  peak  in  Wnet  at  if  =  295  1pm.  Similarly, 
Fig.  7b  shows  the  peak  in  Wnet  with  increases  in  Va.  Increases  in  Va  more  significantly  impact  WyE 
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(c)  Ptube  (cm) 


Fig.  6.  Changes  in  power/cost  ratio  and  net  power  output  with  variation  in  single  parameter  from  optimal  point:  (a) 
/te,  (b)  dmaj  and  (c)  ntube. 


than  increases  in  Vf.  However,  the  dramatic  increase  in  fVTE  with  Va  comes  with  a  rapid  increase  in 
Wfan,  and  Wnel  has  a  sharp  peak  at  the  optimal  point  where  Va  =  2.12  m3/s.  For  all  the  TE  heat 
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(a) 


Vf  (Ipm) 
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(b)  Va  (m3/s) 

Fig.  7.  Changes  in  TE  and  net  power  output  with  variation  in  flow  rate  for  optimized  remote  heat  exchanger:  (a)  V( 
with  associated  pump  power,  (b)  varying  Vn  with  associated  fan  power. 

exchanger  designs,  increases  in  Va  or  V{  increase  the  fan  or  pump  power,  respectively  and,  at  some 
point,  outweigh  the  gains  in  TE  power  generated,  and  thus,  net  power  reaches  a  maximum  at 
intermediate  flows. 


5.  Conclusion 

A  numerical  model  has  been  created  for  an  integrated  thermoelectric  heat  exchanger  in  a  cross 
flow  configuration.  The  features  of  the  model  have  been  validated  with  cross  flow  heat  exchanger 
data  and  with  simple  experiments  integrating  hot  fluid  and  cooling  air  flows  with  Bi2Te3  TE 
modules.  The  results  from  the  validation  studies  confirm  that  the  effectiveness  of  the  system  is 
extremely  sensitive  to  the  operating  temperature  range  of  the  thermoelectric  device.  Using  the 
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validated  model,  a  systems  level  optimization  study  was  conducted  to  show  the  potential  of  using 
cross  flow  heat  exchangers  with  hot  side  fluid  temperatures  characteristic  of  coolant  flows  for  IC 
engines  for  TE  based  waste  heat  recovery. 

The  results  show  that  a  net  power  output  of  1  kW  can  be  achieved  for  a  modestly  sized  heat 
exchanger  core  such  that  the  net  power  density  based  on  heat  exchanger  volume  is  ~45  kW/m3. 
Optimization  for  a  power/cost  ratio  objective  function  with  a  minimum  net  power  requirement  of 
1  kW  indicated  that  the  power  per  cost  can  reach  as  high  as  1.1  kW/$  10,000.  Parametric  studies 
using  the  optimal  solution  as  a  baseline  show  the  importance  of  different  design  variables  for 
waste  heat  recovery  and  further  elucidate  the  important  design  tradeoffs  between  minimizing  the 
air  side  A P  and  optimal  TE  performance.  Further  work  is  needed  to  assess  the  manufacturing 
process  for  integrating  TE  material  into  heat  exchanger  design  in  a  more  cost  effective  manner  to 
reduce  the  current  costs  of  TE  waste  heat  recovery. 
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